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We analyse the results of recent measurements of nonthermal emission from individ- 
ual supernova remnants (SNRs) and their correspondence to the nonlinear kinetic 



Q^i theory of cosmic ray (CR) acceleration in SNRs. It is shown that the theory fits 



these data in a satisfactory way and provides the strong evidences for the efficient 



|H ' CR production in SNRs accompanied by significant magnetic field amplification. 

^ . Magnetic field amplification leads to considerable increase of CR maximum energy 

so that the spectrum of CRs accelerated in SNRs is consistent with the requirements 
for the formation of Galactic CR spectrum up to the energy ~ 10^'' eV. 
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1 Introduction 



The main reason why supernova remnants (SNRs) are considered as a cosmic 
ray (CR) source is a simple argument about the energy required to sustain 
the Galactic cosmic ray (GCR) population against loss by escape, nuclear 
interactions and ionization energy loss. The mechanical energy input to the 
Galaxy from each supernova (SN) is about 10^^ erg so that with a rate of 
about one every 30 years the total mechanical power input from supernovae 
is of the order 10^^ erg/s (e.g. Berezinsky et al. . ll990l ). Thus supernovae have 



enough power to drive the GCR acceleration if there exists a mechanism for 
channeling about 10% of the mechanical energy into relativistic particles. The 
high velocity ejecta produced in the supernova explosion interacts with the 
ambient medium to produce a system of strong shocks. The shocks in turn 
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can pick up a few particles from the plasma flowing into the shock fronts and 
accelerate them to high energies. 



The only theory of particle acceleration which at present is sufficiently well 
developed and specific to allow quantitative model calculations is diffu sive ac- 
celeration applied to the strong shocks associated with SNRs (e. g. see iDrurvl. 



1983 



1991 



Blandford and Eichlerl. 119871 : iBerezhko and Krymskyl . Il988l : iJones and EUisod . 



Malkov and Drurvl . 12001 



for review). Considerable efforts have been 
made during the last years to empirically confirm the theoretical expectation 
that the main part of GCRs indeed originates in SNRs. Theoretically progress 
in the solution of this problem has been due to the development of the ki- 



netic nonlinear the o ry of diffusiv e shock acceleration (IBerezhko et al.l . Il996 



Berezhko and Volkl . 119971 . l2000al ). The theory includes all the most relevant 



physical factors, essential for SNR evolution and CR acceleration, and it is able 
to make quantitative predictions of the expected properties of CRs produced 
in SNRs and their nonthermal radiation. 



Detail information about high-energy CR population in SNRs can be obtained 
from the observations of the nonthermal emission produced by accelerated CRs 
in SNRs. The electron CR component is evident in a wide wave length range 
by the radiation that they produce in SNRs, from radio to 7-ray emission, 
whereas in the case of the nuclear CR component a 7-ray detection is the only 
possibility to find it. If this nuclear component is strongly enhanced inside 
SNRs then through inelastic nuclear collisions, leading to pion production 
and subsequent decay, 7-rays will be produced at the detectable level. 



High-energy 7-rays can also be produced by CR electrons due to the inverse 
Compton (IC) scattering of the background photons. Therefore it is not so 
obvious which CR component (nucl ear or electron) produ ces TeV emission. 



detected from Cassiop e ia A (Cas A) flAharonian et al.l. 120011). RXJ171 3 7-3946 



flMurashi et all. I2OOOI: [Enomoto et al.l . I2OO2I : lAharonian etaD . boOSl . l2006a[ ) 



and Vela Jr ( Aharonian et al. . 2006bl ). For this determination the SNR mag- 
netic field value plays the crucial role. The relative role of electrons in 7-ray 
production is low if SNR magnetic field is strongly amplified. In addition, if 
the magnetic fields in SNRs are considerably larger than the typical interstel- 
lar value, then SNRs are able to generate a power-law spectrum of acceler- 



ated CRs at least up to the knee energy (IBerezhko and Ksenofontovl . 11999 



Berezhko and Volkl . 120041 ). In fact, there is strong evidence that the actual 
magnetic field in SNRs is indeed intrinsically enhanced: the analysis of the 
nonth ermal emiss i on of young SNRs show that all of them have amplified 



fields ( IVolk et al.l . l2005l : iParizot et al.l . l2006l ) . Required strength of the mag- 



netic field, that is significantly higher than the interstellar medium (ISM) 
field, has to be attributed to nonlinear field amplification at the SN shock 
by CR acceleration itself, i n fact by the dominant nuclear CR component 
(iLucek and Bell I2OOOI : Isdl . l2004f l. 
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The application of kinetic theory to individual SNRs (see iBerezhkd . l2005l . for a 
review) has demonstrated its power in explaining the observed SNR properties 
and in predicting new effects like the extent of magnetic field amplification, 
leading to the concentration of the highest-energy electrons in a very thin 
shell just behind the shock. The theory also predicted the enhancement of 
the high-energy galactic 7-ray background radiation due to the contribution 
of CRs confined in parent SNRs, that was partly confirmed by the measured 
background TeV-emission. 

In order to perform the detail comparison of theoretical expectation with the 
experiment one needs sufficient number of individual SNRs with known values 
of relevant physical parameters such as the age, distance, ISM density. Unfor- 
tunately the number of such SNRs are very limited, especially thous which are 
seen in all wavelength from radio to 7-ray . It is why every new experimentally 
established property of nonthermal emission of SNR represents considerable 
interest for theoretical analysis. Therefore we start this review with a brief 
consideration of three SNRs for which essentially new information about the 
efficiency of CR production was recently obtained. Then we shall consider the 
correspondence of the interior magnetic field values extracted from the analy- 
sis of the experimental data with theoretical expectation. Finally, we discuss 
the maximum energy of CRs which can be achieved during their acceleration 
in SNRs, taking into account realistic values of magnetic field in SNRs. 



2 Evidences for efRcient CR production in SNRs 



Experimental and theoretical study of young galactic SNRs, Tycho's SNR, 
RX J1713. 7-3946 and SN 1987A, performed recently provides very useful in- 
formation about CR acceleration by SN shocks. Therefore below we briefly 
analyse these results. 



2.1 Tycho's supernova remnant 



The kinetic theory has been applied in detail to Tycho's SNR, the result of 
a type la SN explosion in a (roughly) uniform interstellar medium (ISM), in 
order to compare the results with existing data, using a stellar ejecta mass 
Mgj = 1.4M r7i, distance d = 2. 3 kpc, and ISM number density Nh = 0.5 H- 



atoms cm ( iVolk et al.l . |2002| . |2005| ). For these parameters a total hydrody- 



namic explosion energy Egn = 0.27 x 10^^ erg was derived to fit the observed 
size Rs and expansion speed Vg. A rather high downstream magnetic field 
strength Bd ~ 300 /xG is nee ded to reproduce the o bserv ed steep and con- 
cave radio spectrum (see also [Reynolds and EUisonl . Il992l ) and to ensure a 
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smooth cutoff of the synchrotron emission in the X-ray region. Overall, very 
good consistency of the predictions of the nonlinear theory with the existing 
observational data was achieved. 



Using Chandra X-ray observations IWarren et al.l (120051 ) have recently esti- 



mated the ratio between the radius Rc of the contact discontinuity (CD), 
separating the swept-up ISM and the ejecta material, and the radius Rs of 
the forward shock. The large mean value Rc/Rs = 0.93 of this ratio was in- 
terpreted as evidence for efficient CR acceleration, which makes the medium 
between those two discontinuities more compressible. Th is result was ana l- 
ysed within already performed kinetic nonlinear approach (j Volk et al.l . l2006l ) . 
Note, that the increase of the ratio Rc/Rs due to the SN shock modification 
by efficiently accelerated CRs was predicted for the case of Kepler's SNR by 



Decourchelle et all (12000f ). 



Fig.l and Fig.2 show the calculations of s hock and CD re l ated quantities 



which were part of the earlier considerations (IVolk et al.l . 120021 . |2005| ). together 
with the azimuthally averaged experimental data available at the time. The 
calculated shock as well as CD radii and speeds are shown as a function of time 
for the two different cases of interior magnetic field strengths = 240 fiG 
and Bd = 360 required by the value of the synchrotron frequency, where 
the spectrum hardens relative to the low-energy radio spectrum. 



According to Fig. la Tycho's SNR is nearing the adiabatic phase. The softness 
of the observed low-energy radio spectrum - relative to a test particle spectrum 
- required a proton injection rate rj = 3 x 10~^. This implies a significant 
nonlinear modification of the shock at the current age of t = 428 yrs. A 
larger magnetic field lowers the Alfvenic Mach number and therefore leads to 
a decrease of the shock compression ratio, as seen in Fig. lb. The result is a 
total compression ratio cr = 5.7 and a subshock compression ratio ag = 3.5 for 
Bd = 240 fiG. In turn a = 5.2, ct, = 3.6, for B^ = 360 fiG. 

Therefore, including CR acceleration at the outer blast wave, the calculated 
value of the ratio Rc/Rs for Sd = 360 /iG is slightly lower than for 5^ = 
240 /iC. At the current epoch we have Rc/Rs ~ 0.90 which is somewhat lower 
than the value Rc/Rs = 0.93 inferred from the Chandra observations. 



Projecting a highly structured shell onto the plane of the sky tends to favor 
protruding parts of the shell. Therefore the average radius measured in pro- 
jection is an overestimate of the true average radius. A nalyzing the amount 



of bias from the projection for the shock and CD radii IWarren et al.l (120051 ) 



found a corrected "true" value Rc/Rs = 0.93 which is lower than their mea- 
sured "projected average" value Rc/Rs = 0.96. One needs to make such a 
purely geometrical correction independently of the nature of the factors which 
produce the shell structures, if one measures the size of the quasispherical 
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structured shell in projection instead of dealing with a 2-dimensional shell 
crossection. In the latter case there would be no need for such a correction. 



Irrespective of this need to correct the observations of the CD position, if one 
starts from a spherically symmetric calculation of the CD radius as we do, one 
has to take into accou nt that the actual CD is subject to the Ray leigh- T aylor 
(R-T) instability (e.g. IChevalier et all Il992l : IWang and Chevalierl . I2001I ). and 
thus a correction is needed to compare such a 1-D calculation with observations 
of the CD. In the nonlinear regime the instability leads to effective mixing of 
the ejecta and swept-up ISM material with "fingers" of the ejecta on top 
of this mixing region. The latter has a radius larger than Rj,^. The minimum 
correction to Rl^ is the radial extent ARc of the mixing region above Rl^. We 
estimate ARr ~ 0.51, where / ~ 0.1 -R^^ is the longest finger size according to 



the numerical modeling of IWang and Chevalied (120011 ) (albeit without particle 
acceleration). This leads to a rough estimate of the corrected CD radius R'^ = 
l.05R}^. 



As it is seen in Fig. 2, the c omparison o f the correct ed values R'^/ Rs, according 



2OO2I . I2OO5I ). with this experimentally 



to the earlier calculations fiVolk et al 
estimated value Rc/Rs = 0.93 shows quite good agreement. 



The agreement between theoretical solutions for effective particle accelera- 
tion with the measurements of the shock and discontinuity radii has to be 
considered as new evidence for strong nonthermal effects in Tycho's SNR. 
New Northern Hemisphere TeV detectors is expected to detect this source at 
TeV-energies in, predominantly, hadronic 7-rays at an energy flux level above 
2 X 10~^^erg cm~^ s~^. As a corollary the detection of a TeV signal is not only 
important by itself, but it is also crucial for the correct determination of all 
other key SNR parameters. 



2.2 Supernova remnant RX J1713-3946 



RX J1713-3946 is a shell-type supernova remnant (S NR), located in the Galac- 
tic pl ane, that was discovered in X-rays with ROSAT (jPfeffermann and Aschenbach 
19961) . Further stud y of this SNR wit h the ASCA sateUite by iKoyama et al 



( 119971 ) and later by ISlane et al.l (119991 ) have shown that the observable X-ray 
emission is entirely non-ther mal, and this propert y was confirmed in later 



XMM-Newton observations (ICassam-Chenai et al.l . 120041 ). The radio emis- 



sion of this SNR is weak: only part of the shell could be detected in ra- 
dio synchrotron emis sion up to now, with a poorly known spectral form 
( jlazendic et all . 12004 ). 



RX J1 713-3946 was a l so detected in very high ene rgy 7-rays wit h the CANGA- 
ROO (iMurashi et al.l . l200d[Enomoto et al.l . l2002h and H.E.S.S. (lAharonian et al. 
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2006al ) telescopes. Especially the latter observations show a clear shell struc- 
ture at TeV energies which correlates well with the ASCA contours. 



The difficulty for the theoretical description is the fact that several key param- 
eters of this source are either not known or poorly constrained. This already 
concerns the distance and age of the obje ct. It was demonstrated th at consis- 
tent description of this object is achieved (IBerezhko and Volkl . l2006l ) following 
present consensus which puts the distance at 1 kpc, the age to about 1600 years 
and that the primary explosion must have been a type Il/Ib SN event with a 
massive progenitor star whose mass loss in the main sequence phase created 
a hot wind bubble in a high-density environment. The solution for the overall 
remnant dynamics then yields the value for the expansion velocity of the outer 
shock, given the total mechanical energy E^^ = 1.8 x 10^^ erg released in the 
explosion. To obtain a consistent solution for the broadband nonthermal emis- 
sion the injection rate r] = 3 x 10^^ and interior magnetic field ~ 130 /xG 
are needed. 



The calculated overall broadband spectral energy distribution is displayed in 
Fig. 3, together with the experiment al data from ATCA at ra dio wavelengths, 
as estimated for the full remnant by lAharonian et al.l (l2006al ). the X-ray data 
from ASC A, and the most r ecent ly measured TeV gamma-ray spectrum from 
H.E.S.S. (lAharonian et al.l . 120071 ). The overall fit is impressive, noting that 
the choice of a few key parameters like rj, B^, -Esn in the theory allows a spec- 
trum determination over more than 19 decades. Note, that calculated gamma- 
ray stectrum (Berezhko and Volk, 2006) agrees with the new H.E.S.S. data 
(Aharonian et al., 2007) even better than with earlier H.E.S.S. measurements 
(Aharonian et al., 2006a). 



The prop erties of small scale st ructures of SNR RX J1713. 7-3946 s een i n 
X-rays by lUchiyama et al.l (120031 ) . and in particular by iHiraga et al.l (120051 ) . 



provide even stronger evidence that the magnetic field inside the SNR is indeed 
considerably amplified. 



In order to find out if the filamentary structure found by (IHiraga et al.l . |2005| ) 
is indeed consistent with the efficient CR production projected radial proffie 
J(ej,, p) calculated for the X-ray energy = 1 keV is shown in Fig. 4. The 
theory predicts the peak of the emission just behind the shock front with 
projected thickness Ap/Rg ^ 10~^ that corresponds to an angular width Aip ^ 
0.4'. This width is significantly thinner than the observed width of Aip ^ 2.2'. 
However, one should take into account the smoothing procedure with which 
the data were obtained/presented. If we apply a point spread function of width 
= 12.8", corresponding to the smoothed XMM data, to the calculated 
proffie, we obtain the broadened proffie , which is shown in Fig.4 by the dashed 
curve. It has a width Ap/ R^ ~ 6.6 x 10~^, or A ijj f» 2', and is consistent with 
the observational proffie of IHiraga et al. (l2005h . 
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In addition the projected radial profile (dash-dotted line) which corresponds to 
the test particle limit is presented in Fig.4. In this case the proton injection rate 
was taken artificially so small [r] = 10~^), that accelerated CRs do not produce 
any significant nonlinear effect, like shock modification and/or magnetic field 
amplification. Therefore interior magnetic field value = 20 fiG is much 
lower in this case. Due to this fact synchrotron losses are not relevant on 
account of the low magnetic field, the peak of the emission is broader by a 
factor of five, inconsistent with the XMM observations. 

The projected radial 7-ray emission profile, calculated for the energy = 
1 TeV, is presented in Fig. 5. As a result of the large radial gradient of the gas 
density and the CR distribution inside the SNR, the theoretically predicted 
radial profile of the TeV-emission is concentrated within a thin shell of width 
Ap ^ 0.1i?s- Since the H.E.S.S. instrument has a finite angular resolution we 
present in Fig.4 also the modified radial profile smoothed with the Gaussian 
point spread function with dp = Ap = 0.08-Rs, that corresponds to the angular 
resolution 2cr^ = 0.1°. 

As shown by Fig. 5 the smoothed radial profile of the TeV-emission is much 
broader and - what is most interesting - it is characterised by a maximum to 
minimum intensity rati o J™'" /J™'" = 2. 2 . Such a ratio is consistent with the 



H.E.S.S. measurement (lAharonian et al.l . l2006al ) which obviously gives only a 



lower limit to the sharpness of the 7-ray profile. Note that the broad radial 
profile of the TeV-emission with J^'^^/J™™ ^ 2 measured by the H.E.S.S. 
instrument is indirect evidence that the actual radial profile is significantly 
sharper, with a higher ratio J™'*''/J™" > 2.2. 

We conclude that the present observational knowledge of SNR RX J1713.7- 
3946 can be interpreted by a source which ultimately converts about 10% 
of the mechanical explosion energy into nuclear CRs and that the observed 
high energ y 7-ray emission of SNR RX J1713. 7-3946 is of hadronic origin. 



Note, that iMalkov et al.l (120051 ) have obtained the same type of conclusion 



sugg esting a rather different 7 -ray spectrum formation scenario, compared 



with iBerezhko and Volkl ((20061). 



2.3 Supernova remnant SN 1987 A 



Supernova 1987A occurred in the Large Magellanic Cloud. It has been exten- 
sively studied in all wav elengths frorn radio to 7-ray . The initial short out- 
burst of radio emission (ITurtle et al.l . 119871 ) is attributed to the synchrotron 
emission of electrons accelerated by the SN shock propaga ted in the free wind 



of pre supernova star, which was the blue supergiant (BSC) (I Chevalier and Fransson 



19871 ) . After about 3 years radio emission was detected again (IStaveley-Smith et al. 



7 



1992 



Gaensler et all ll997D a s well as the monotonic ally increased X-ray emis- 



sion (IGorenstein et al.l . Il994l : iHasinger et al.l . Il996l ). This second increase of 



emission is attributed to the entrance of the outer SN shock into the ther- 
malized BSG wind and then in the H II region occupied by much more dense 
matter, consists of the swept up w ind of red supergiant (RSG) progenitor star 



( iChevalier and Dwarkadad . 119951 ) 



The advantage of SN 1987A for theoretical consideration is that the key pa- 

lOMp 



05 



rameter values are very well determined: stellar ejecta mass M, 
dista nce d = 50 kp c, hydrodynamic explosion energy Egn = 1-5 x 10^^ erg 
(e.g. iMcCrayl . 119931 ). During an initial period the shell material has a broad 
distribution in velocity v. The fastest part of this ejecta distribution is de- 
scribed by a power law dMej/dv oc with k = 8.6. 

The kinetic nonlinear model for CR acceleration in SNRs has been in de- 
tail applied to SN 1987A, in order to compa re its results with observational 
properties (IBerezhko and Ksenofontovl . l2006l ). 



A rather high downstream magnetic field stren gth > 1 mG is needed 



to re produce the observed steep radio spectrum (IBerezhko and Ksenofontovl . 
20001 ). The required strength of the magnetic field have to be attributed to 
nonlinear field amplification at the SN shock by CR acceleration itself. Since 
for all the thoroughly studied young SNRs, the ratio of magnetic field energy 
density B^/Sn in the upstream region of the shock p r ecurs or to the CR pres- 
sure Pc is about i?Q/(87rPc) ~ 5 x 10~^ (IVolk et al.1 . l2005l ) and CR pressure 
in young SNRs has a typical value Pc ~ 0.5^0^^, the upstream magnetic field 
Bq = Bd/cr was taken in the form 



Bn 



2n X 



Calculated shock radius Rg and speed Vg shown in Fig. 6a as a function of 
time are in satisfactory agreement with the values obtained on the basis of 
radio and X-ray measurements. Note that radio data compared with X-ray 
data gives larger shock size at any given time t and calculated radius Rs{t) 
goes between these two sets of data. 

To fit the spectral shape of the observed radio emission a proton injection 
rate rj = 3 x 10~^ is needed. This leads to a significant nonlinear modification 
of the shock: as it is seen in Fig.6b total shock compression ratio a ~ 5.3 
is essentially larger and a subshock compression ratio cr^ ~ 2.8 is lower than 
classical value 4. 

Since during the last 15 yrs the shock speed and the gas density are almost 
constant CR energy content grows roughly linearly with time (see Fig. 6c). 
It gives the natural explanation of the linear increase of the radio emission 
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detected during this evolutionary period. 



Strongly modified SN shock generates CR spectrum oc , which is very 
steep at momenta p < rUpC, with index 7 = {as + 2) /{as — 1) ^ 2.7. CR 
electrons with such a spectrum produces synchrotron radioemission spectrum 
S,y OC z/^" with spectral index a = (7 — l)/2 ^ 0.9, that very well corresponds 
to the experiment, as it is seen in Fig. 7, where synchrotron energy spectra 
uSy, calculated for five subsequent epoch together with the experimental data 
are presented. Note that CR spectrum has a concave shape: it becomes flatter 
at higher momenta p. As a consequence synchrotron spectrum S^{i') is also 
concave as it is clearly seen in Fig. 7 a.t u < 10^^ Hz. Radio data reveal this 
property in good consistency with theoretical prediction. 



Strong downstream magnetic field ^ 15 mC, that corresponds to the up- 
stream field -Bo ~ 3 mG (see Fig. 7), provides synchrotron cooling of electrons 
with momenta p > lOrripC that in turn makes synchrotron spectrum at high 
frequencies u > 10^^ Hz very steep (see Fig. 7). Concave shape of electrons con- 
tinuously produced at the shock front together with their synchrotron cooling 
lead to a formation of two peaks in synchrotron energy spectrum uS^,. The 
first one at z/ ^ 10^^ Hz corresponds to CR electron momentum p ^ lOrripC 
above which synchrotron energy looses are relevant, whereas the second peak 
at 1/ ~ 10^^ Hz corresponds to the maximum momentum p ^ W^rripC of 
accelerated electrons. Under this condition calculated synchrotron flux at fre- 
quency u ~ 10^'' Hz, which corresponds to the photon energy = 0.5 keV, 
is below the measured flux at the epochs t > 3000 d. Since the contribution 
of the nonthermal radiati on in the obse r ved X -ray emission of SN 1987A is 



not very well known (e.g. iMichael et al.l . |2002| ). the observed X-ray flux has 



to be considered as the upper limit for the expected nonthermal emission. At 
early epoch t < 2500 d however the calculated flux exceeds the measured one 
(see the curve, corresponding t = 1970 d in Fig. 7). This can be considered as 
indication, that the actual magnetic field Bq is few times larger then given by 
the Eq.(l). 



At lower magnetic field B^ ^ 2 mG synchrotron losses of high energy CR 
electrons are considerably smaller compared with the previous case. Due to 
this fact synchrotron s pectra considerably exceeds at an y given epoch the 
measured Chandra flux (IBerezhko and Ksenofontovl . |2006| ) . Therefore one can 
conclude, that the actual interior magnetic field strength is not lower than 
5 mC. 



We see that quite a reasonable consistency of kinetic nonlinear theory with 
most of the observational data of SN 1987A is achieved, that gives the evidence 
of the efficient CR production leading to a strong shock modification and 
strong magnetic field amplification. 
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Calculated 7-ray integral energy flux at all energie s is dominated by the vr"- 
decay component (IBerezhko and Ksenofontovl . l2006l ) . At the current epoch the 
expected 7-ray energy flux at TeV-energies is about e^F^ ^ 2x 10~^^ erg/ (cm^s) 
and during the next four years it is expected to grow by a factor of two. The 
detection of high energy 7-rays from SN 1987A, which could be done by in- 
struments like H.E.S.S., would provide very clear evidence for the efficient 
acceleration of nuclear CRs. 



2.4 Discussion 



Detailed consideration performed within a frame of nonlinear kinetic theory 
demonstrates that the CR production efficiency in young SNRs is consistent 
with requirements for the Galactic CR sources. Recent experimental and theo- 
retical studies of Tycho's SNR, RX J1713. 7-3946 and SN 1987A provides new 
very strong evidence for that. 



Consideration performed by iHillad (120051 ) led him to the opposite conclusion: 
that TeV 7-ray emission from the mostly fairly young SNRs (SN 1006, Cas A 
and Tycho's) is unexpectedly low. Note, however, that this consideration is 
based on the simplified estimates, which does not take into account such an im- 
portant aspect as the necessity of renormalization of CR production efficiency, 
predicted in spherical approach. We have argued before that ion injection is 
quite efficient at the quasiparallel portions of the shock surface and it is ex- 
pect ed to be strongly suppressed at the quasiperpendicular part of the shock 



I see IVolk et al.l . |2003| . for details). Therefore one should renormalize the re- 
sults for the nucleonic spectrum, calculated within the spherically symmetric 
model. The lack of symmetry in the actual SNR can be approximately taken 
into account by a renormalization factor /re < 1, roughly /re = 0.15 to 0.25, 
which diminishes the nucleonic CR production efficiency, calculated in the 
spherical model, and all effects associated with it. 



For SN 1006 which appears to approximate this ideal case, efficient CR pro- 
duction is expected to arise within two polar regions, where the SN shock 
is quasi-parallel. The calculated size of the efficient CR production regions 
which amounts to about 20% of the shock surface corresponds ver y well to the 



observed sizes of the bright X-ray synchrotron emission regions (jVolk et al. 



20031 ). This implies a renormalization factor of 0.2. In addition, the distance to 
the object d = 1.1 kpc used by iHillasI (2005) for the case of SN 1006 is consid- 
erably lower than the most reliab le determination of the distance d = 2.2 kpc 
to SN 1006 (IWinkler et al.l . l2003h . 



Due to t he above two factors the estimate of TeV 7-ray flux from SN 1006, 
made by iHillad (120051 ) . should be reduced by a factor of 20, so that it will 
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be in agreement with H.E.S.S. upper limit. Note, that the study of SN 1006 
on the basis of kinetic nonhnear theory demonstrated that the lack of a TeV 
signal that follows from the non-detection by the H.E.S.S. instrument does 
not invalidate the theoretical picture which gives a consistent description of 
the nonthermal emission characterist ics; it rather implie s a co nstraint on the 



ambient gas density A^h < .1 cm ^ (IKsenofontov et al.l . |2005| ) . Using XMM- 



Newton X-ray observations lAcero and Balletl (120061 ) recently confirmed that 
the ISM density is indeed low, iVn < 0.1 cm~^. The example of SN 1006 
illustrates that real difficulty in reliable prediction of expected 7-ray emission 
from SNRs is not a theoretical one, but rather due to poorly known values of 
relevant astronomical parameters, like di stance to the object and a mbient ISM 
density. /^From theoretical consideration (IKsenofontov et al.l . |2005| ) it becomes 
clear that in the case of SN 1006 and in other similar cases 7-ray measurements 
will at the same time give a reliable estimate of the ambient ISM density. 



Similar sit uation i s for the two other cases of Cas A and Tycho's SNR, con- 
sidered by iHillad (120051 ): renormalization of his estimates, made in spherical 
approach, and the adoption the most reliable astronomical parameters provide 
the value of the expec ted 7-ray f lux co nsistent with the existing measurements 



(or upper limits) (see iBerezhkd . l2005l ). 



3 Amplified magnetic field in SNRs 



All young SNRs (including three above cases) to which nonlinear kinetic ac- 
celeration theory has been successfully applied up to now exhibit the strongly 
amplified niagnetic field as result of very efficient acceleration of nuclear CRs 



( jVolk et al.l . l2005l : lBerezhkd . I2OO5I ). Magnetic field extracted from the spectral 



properties of the overall synchrotron emission and from the observed filamen- 
tary structure of the nonthermal X-ray emission agree very well and their 
value can be expressed in the form 



5oV(87r) = 5 X lO^^Pe 



(2) 



where Pc is CR pressure at the shock front, Bq ^ Bjsm is the upstream field, 
which is already amplified and therefore considerably exceeds magnetic field 
value BisM ~ 5 fiG in the ambient ISM. Let us consider the consistency of 
such a field strength with the plasma physics. 



Several attempts have been undertaken to model highly nonlinear turbu- 
l ent magnetic field, created by accelerated CR s upstream of the s trong shock 



( Bell and Lucekl . l2001l : IVladimirov et al.l . l2006l : lAmato and Blasil . 120061 ) . Two 



aspects of this process are of significant importance: the amount of magnetic 
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field amplification and the extent of of gas heating due to the turbulence dis- 
sipation. Unfortunately up to now it is not clear in which proportions CR 
working on the gas in the precursor region is divided between gas heating and 
magnetic field amplification. Therefore we restrict ourselves here by a simple 
estimates in order to demonstrate that considerable magnetic field amplifica- 
tion Bq ^ BisM, consistent with Eq.(2) is indeed possible upstream of the 
strong shock. 



There are two different mechanisms operating upstream of the fast super- 
Alfvenic shock and leading to the amplification of the chaotic magnetic field. 
The first one is related with the resonant interaction of CRs with the Alfven 
waves, which due to the CR anis otrop y trigge rs an instability, fol l owed by the 
intense Alfven wave excitation ( Belli . I1978| : ,McKenzie and Volkl . Il982l ). The 



other one, recently proposed by iBelll (120041 ). is non- resonant amplification of 
MHD perturbations due to Lorentz force, created by CR current. 



Integral effect of resonant A lfven wave excitation can be written in the form 
flMcKenzie and Volkl . [l982h 



2m- 



B 



dE_ 
dx 



dPc 
dx 



(3) 



if we suggest that the CR working CAdPc/dx goes on the excitation of the 
Alfven waves. Here Eb = (5-B)^/(87r) is magnetic component of the Alfven 
wave energy density, ca = B/y/Airp is Alfvenic speed, u is the plasma flow 
speed directed along the x axis to the shock situated at the point x = 0. For 
the purpose of estimates we suggest quasistationary condition and neglect the 
gas speed variation u{x) in the upstream region x < due to the shock mod- 
ification by CR backreaction, taking u = Vg. Assuming that the far upstream 
field Bq = BjsM coincides with the ISM field Bjsm and performing the in- 
tegration of the above equation we have magnetic field energy density at the 
shock front 

Eb/Ebo = MAoPJipu^), (4) 



where Mao = Vs/cao is Alfvenic shock Mach number, Ebo = Bq/^Sh) is 
the energy density of the far upstream field. Since the strong shock provides 
efficient CR acceleration with Pc/{pu^) ~ 1 the above expression shows that 
Alfven wave amplitude grows to a large value 5B ^ Bq in the case of strong 
shock with Mach number Mao ^ 1- 

Resonant magnetic field amplification should be supplemented by two effects: 
non-resonant instability and the Alfven speed variation ca{B) across the pre- 
cursor due to the increase of the mean effective field B. In this case (see 
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Pelletier et al.l . l2006l . for details) instead of (4) we have 

Eb/Eb^ = MAoiB,,r/Bo)Pc/ipu'), (5) 

where 



Bnr=Bo,l^^^ (6) 



is the mean field amplified non-resonantly. Here $ = \og{pmax / f^pc) , Pmax is 
the maximum CR momentum. 

Substituting $ = 10 in the expression (5) we have 

Eb = Sx 10-2(^/3000 km/sy/^Pc. (7) 



Since for all known young galactic SNRs the shock speed is about K ~ lO^^c, 
this expression shows, that existing mechanisms is able to amplify magnetic 
field within the precursor to the level, which considerably exceeds the required 
value (2). Therefore we can suggest that the main amount of the energy given 
by the expression (5) in fact goes on the gas heating due to dissipation of 
MHD perturbations and only small fraction of it, given by (2), indeed goes on 
the magnetic field ampli fication. Numerical simulation of considered process 



performed by iBelll (120041 ) confirmed, that the gas thermal energy an order of 



magnitude larger than the energy content of amplified magnetic field. 

Note that the process of magnetic field ampl i fication is not in cluded in kinetic 
nonlinear theory (e.g. iBerezhko et al.l . Il996l : iBerezhkol . liooi ). It simply pos- 



tulates the existence of already amplified far upstream field Bq, which value 
is consistent with experimental requirements, given by Eq.(2). The consis- 
tency of such a picture is provided by the fact, that the spectrum of CRs 
produced by strong modified shock is very hard so that CRs with highest 
energies have a largest contribution in their energy content. These the most 
energetic CRs produce field amplification on their spacial scale that is the 
precursor size. Therefore CRs with lower energies already 'see' the amplified 
field Bq 3> BjsM- Note also, that kinetic nonlinear model contains in the 
explicit form the suggestion, that the whole CR working goes on the gas heat- 
ing presumably due to some kind of dissipation mechanism. Since according 
to the above estimates, the magnetic field energy which can be potentially 
achieved in the precursor region is considerably larger than required energy, 
given by Eq.(2), it can be considered as indirect evidence for the existence of 
such dissipation processes, which efficiently converts the energy of magnetic 
field perturbations into the gas thermal energy. 
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The existence of nonadiabatic gas heating within the precursor plays an im- 
portant role for the shock s tructure and dynamics. It restricts the shock corn- 
pression ratio by the value (IBerezhko et al.l . 1 19961 : iBerezhko and EUisonl . Il999l ) 



a 



1.5M' 



3/8 



AO ? 



if the whole CR working caAPcI goes on the gas heating. For the appropriate 
upstream magnetic field values, consistent with Eq.(2), it provides moderate 
shock modification with a ~ 6. In the opposite case, when within the precursor 
there are no any addition al gas heating except adiabatic one, the expected 



shock compression ratio is (IBerezhko et al.l . Il996l : IBerezhko and EUisonl . Il999l ) 



1.3M: 



3/4 



50 ' 



(9) 



where Msq = Vg/cso is the sonic Mach number, cso is the far upstream sound 
speed. For the typical ISM parameters this expression leads to much larger 
shock compression a ^ 6. 

In the case when CR working in the upstream region goes on the gas heating, 
the expected gas thermal energy just ahead of the subshock is 

Eg ^ PJMao. (10) 



Appropriate magnetic field value given by Eq.(2) leads to a typical shock 
compression ratio a = 6, that according to Eq.(8) corresponds to Maq ~ 40, 
that in turn according to Eq.(lO) provides gas thermal energy Eg 2.5 x 
10~^Pc- This value Eg very well corresponds to what is expected from the 
expression (7), if indeed main amount of magnetic field perturbation energy 
is dissipated. We therefore conclude that kinetic nonlinear theory provides 
the consistent description at least at shock speeds Vs ~ 3000 km/s, that 
correspond to the most active phase of the SNR evolution that is the beginning 
of the Sedov phase. 



4 Maximum energy of CRs 



The most important consequence of magnetic field amplification in SNRs is the 
substantial increase of the maximum energy of CRs emax = PmaxC, accelerated 
by SN shocks. Let us briefly analyse what value of emax can be expected taking 
into account the amplification of SNRs magnetic field. 



According to theoretical co nsideration (Lucek and Bell . 2000l ) and the analysis 



of the observational data (jParizot et al.l . l2006l ) CR diffusion in the amplified 
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magnetic field is close to the Bohm limit. Therefore cutoff momentum of CRs 
accelerated at an y given SNR evolutionary stage can be expressed in the form 
f lBerezhkol . [l996l ) 



Pm = mpcRsVs/[AK{mpc)], 



where = pc^ /{SZcBq) is CR diffusion coefficient, corresponding to Bohm 
limit; Z is the charge number of CR nuclei. The value of parameter A depends 
on the shock expansion law Rs{t) and on the downstream gas flow radial 
profile w{r). Note that the main factors which determines the value of CR 
cutoff momentum pm are: i) finite shock size; ii) shock deceleration and iii) CR 
adiabatic cooling in the expended downstream region r < -R,, . These factors are 



more relevant compared with the time factor, studied by lLagage and Cesarsky 
(119831 ). that in turn prev e nts t he complete shock modification due to CR 



backreaction (jPrury et al.l . Il995l ) 



First consider the most simple case of uniform ISM, that is typical for type 
la SNe. In such the main part of CR spectrum formed during the SNR 

evolution is produced at the Sedov phase. Since the CR cutoff momentum at 
the Sedov phase Pm oc RsVsBq oc RgV^ oc t~'^^^ decreases with time t the most 
energetic CRs in the power law part of this spectrum are produced at the 
very end of the free expansion phase. Therefore the maximum CR momentum 
Pmax = niax{pm(^)} depcuds on physical parameters as Pmax oc Vq-Ro-BoI^o), 
where Vq, Rq and to are mean ejecta speed, sweep up radius and sweep up 
time respectively. Taking into account th e value of upstrean i magn etic field 



(2) and assuming Pc = O.SpoK^, we have (IBerezhko and Volkl . |200J) 



Pn 



rUpC 



3 X lO^Z 



M, 



-2/3 



1051 erg / \ 1.4Mg 



AT, 



H 



0.3 cm 



-3 



1/6 



'12) 



It has a strong dependence on the SN parameters Esn and M^j, but is only 
weakly dependent upon the ISM density po = lANnmp. 

Maximum momentum of protons, given by Eq.(12) for Esn = 10^^ erg, Mgj = 
IAMq and = 0.3 cm~^ is Pmax ~ 3 x lO^mpC. This required condition to 
reproduce the spectrum of GCRs up to the knee energy is fulfilled due to the 
magnetic field amplification. 

For CR iron nuclei Pmax ~ lO^m^c that provides the possibility for the forma- 
tion of GCR spectrum up to the energy ~ 10^^ eV inside SNRs. 

Since during the early free expansion SNR evolutionary phase the SN shock 
speed Vg = {3 — A) X km/s is considerably higher than on the subsequent 
Sedov phase {Vg = (3 — 10) x 10^ km/s), one could expect the production 



15 



of CRs with energies p >_Pmax essentially higher than produced on Sedov 
phase (iBell and Lucekl. 120011) . Let us brie fly consid er this possibility, which 
was s tudied by iBerezhko and Volkl (120041 ) and by iPtuskin and Zirakashvili 
fl2005h . 



During the free expansion phase the SN shock evolution is determined by the 
velocity distribution of the fastest part of the ejecta which is described by a 
power law 



dMej/dv oc V 



2-k 



(13) 



with A; = 7 to 12 fl Jones et al.l . ll98ll : IChevalie^ . ll982h . The ejecta kinetic energ; 



y 

contained in the fraction with speed v within the interval dv is dEej{v) oc 
v'^'^dv. The CR spectrum created by the strongly modifled shock is quite 
hard: N{p) oc p"'*', 7 < 2. For the sake of a rough estimate we assume that 
all the above ejecta energy goes into CRs with maximum momenta Pm(K), 
where the shock speed is roughly the ejecta speed K ~ v. This gives 



N{pm)Pmdpm OC dE, 



(14) 



in the ultrarelativistic case. Taking into account that Pm oc RgVgB oc V^, we 
obtain 

7 = 2 + (A; - 5)/a. (15) 



Since in the case of uniform ISM R, oc t^^-^)/^ we have a = (9 - k)/3 if 
we suggest for the simplicity Bq oc poV^. Even i n the case of type la SNe 
with a typical value k = 7 Eq .(15) gives 7 = 5 (IBerezhko and Volkl . |2004J : 



Ptuskin and Zirakashvili! . |2005| ). It means that the spectrum of CRs with mo- 
menta p > Praax accelerated during the free expansion phase is so steep, that 
it hardly plays any role for the formation of the GCR spectrum. 



In th e case of type Il/Ib SNe, which are more numerous in our Galaxy (ITammann et al 
19941 ). during the initial SNR evolutionary phase SN shock propagates through 
the progenitor star wind. Therefore in this case ambient gas density at the 
shock front varies as po R-'^ that gives a = 2. Even at the lowest value 
= 8 the resulting CR spectrum is very steep, 7 = 3.5. 

We therefore conclude, that CRs accelerated during the free expansion phase 
do not play a role for GCR spectrum, except the case if the actual value of the 
parameter k is as small as < 5.4 for type la SNe and k < Q for type Il/Ib 
SNe or the GSM has more complicated structure than considered above. 



In Fig. 8 we present GCR energy spectra, expected in the Galaxy, assuming 
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that S NRs are the main GCR source. Compared with our previous calcula- 
tions (IBerezhko and Ksenofontovl . Il999l ) much higher magnetic field in SNR 
is used, which varies during the SNR evolution according to Eq.(2). To obtain 
these spectra J(efc) oc J s{^k)^esc{.R) the corresponding source spectra Jsi^k), 
produced in SNR during its evolution, were calculated for a typical set of rele- 
vant parameters {Egn = 10^^ erg, M^j = IAMq, Nh = 0.3 cm~^), taking into 
account their steepening due to the GCR escape from the Galaxy with a time 
scale Tp,^r. oc R~^, where R(ek) is t he rigidity, is the kinetic energy of CR 
particle (IBerezhko and Volkl . 120071 ) . 



For a comparison with the e xperiment we repres ent in Fig. 8 the data of only 



three experiments, ATI C-2 (IPanov et al.l. 120061) . JACEE (ITakahashi et al. 
19981 ) and K AS CADE (lAntoni et al.l . |2005| ). which on our view are repre- 



sentative in the corresponding energy ranges. One can see, that the theory 
in a satisfactory way fits the existing data up to the energy ~ lO"*^"^ eV. 
The e ssential exception is the measured in the recent ATIC-2 balloon experi- 
ment (jPanov et al.l . |2006| ) helium spectrum, which opposite to the theoretical 
expectation is noticeably fiatter than the protons spectrum. Note that the 
spectral shape of protons and heli um measured in the previous ATIC-1 ex- 
periment at energies 10 



10^ GeV flAhn et al 



20061 ) are very similar. Similar 
form of the spectra for different GCR nuclei at energies ek = 10^ — 10® GeV 
if observed could be considered as the indirect confirmation of GCR origin 
in SNRs. Therefore the confirmation of ATIC results for protons and helium 
by other experiments with comparable or even better precision is very much 
needed. 



Another difficulty for the theory is the fact that to get a consistency with the 
observed GCR spectra the value = 0.75 is neede d, which is a bit beyond the 



experimentally determined interval fi = 0.3 — 0.7 ( Berezinsky et al. , 1990l ) 



In addition, strong energy dependence of the GCR escape time Tesd^) makes it 
difficult to interpret the observed CR anis otropy , whic h is very low, less than 
one percent at energies e = 10^^ — 10^^ eV (jHillasi . l2005l ). Here e is total energy 
of CR particle. Total fiux of GCRs escaping from the Galaxy oc e~'^ has 
energy dependence hke CR source spectrum Jg oc e^'^ with •y ^ 2. Therefore 
GCR anisotropy outside the region where CR sources are situated ( Galactic 
disk volume) increases with energy Fc/J oc e^, like inverse residence time 
T~l. If fi is as big as 0.6 the expected anisotropy at e = 10^^ eV exceeds 



the level of 100% (jHillasi . l2005l ). Since our observing positi on is not out side 
the source layer, we should see a smaller GCR anisotropy (jHillasi . l2005l ). In 
addition one can also suggest that in the region, occupied by CR sources, GCR 
diffusion coefficient is higher than in outer part of GCR confinement volume, 
due to enhanced level of turbulence produced by GCRs streaming away from 
the Galaxy. In such the GCR anisotropy in the source layer could be 

sufficiently small and consistent with the observations. 
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Note that the theory predicts concave spectra at energies e > 10 Z GeV for 
individual CR species and there are some experimental evidences that the 
actual GCR spectra are indeed concave, eventhough one needs more precise 
measurements at e > 10'^ GeV in order to make more strict conclusion. At the 
same time calculated all particle spectrum has almost pure power law form 
up to the knee energy (see Fig. 8). 



According to Fig.8 the knee in the observed all particles GCR spectrum has 
to be attributed to the maximum energy of protons, produced in SNRs. The 
steepening of the all particle GCR spectrum above the knee energy 3 x 10^^ eV 
is a result of progressive depression of contribution of light CR nuclei with 
the energy increase. Such a scenario is confirmed by KASCADE experiment 
which shows re latively sharp cutof f spectra of GCR species at energies e^ax ~ 
3Z X 10^^ eV dAntoni et all . l2005h . At energy e ~ 10^^ eV GCR spectrum is 
dominated by iron nuclei contribution. 



According to the analysis, performed by lAloisio et al.l (|2006( ) the contribution 
of Galactic CR sources into the observed GCR spectrum dominates up to 
the energy 10^'' eV above which it progressively goes down so that at energy 
10^^ eV the contribution of extragalactic sources becomes dominant. At so- 
called second knee energy = 5 x 10^'' eV these two sources provide roughly 
equal contributions into GCR spectrum. As it is seen from Fig.8, calculations 
provide a good fit of the observed GCR up to the required energy = 10^^ eV, 
that is a strong argument for SNRs to be considered as the main Galactic CR 
source. 



Essentially different scenario of the formation of GCR spectrum and the first 
knee in particular was proposed and discussed b y Erl ykin and Wolfandale in 
a series of papers (see Erlykin and Wolfendale . 200Cll . and references there). 



Based on the analysis of measurements performed by ground based installa- 
tions they claim the existence of complex feature of GCR spectrum in the knee 
energy region. This leads them to a conclusion that observed GCR spectrum 
consists of smooth background spectrum (without the knee), produced by 
many CR sources, and superimposed single source (presumably nearby SNR) 
contribution, that forms a bump on the smooth background spectrum usually 
interpreted as a (first) knee. The contribution of the single nearby SNR, lead- 
ing to some peculiarity of GCR spectrum at the detectable level in the knee 



region, is indeed quite possible (IBerezhko and Ksenofontovl . Il999l ). However, 
the existing measurements of CR composition demonstrate, that the knee is 
indeed peculiar point in GCR spectrum, where not only the shape of GCR 
spectrum but also their chemical composition undergoes a noticeable change, 
that is consistent with CR production in typical SNRs, as described above. 
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5 Summary 



Detailed consideration performed within a frame of nonlinear kinetic theory 
demonstrates that the CR production efficiency in young SNRs is consistent 
with requirements for the Galactic CR sources. Recent experimental and theo- 
retical studies of Tycho's SNR, RX J1713. 7-3946 and SN 1987A provides new 
very strong evidence for that. 



Synchrotron radiation from young SNRs provides evidence that efficient CR 
acceleration followed by the significant shock modification and magnetic field 
amplification takes place there. Magnetic field amplification leads to the con- 
siderable increase of maximum energy of CRs, accelerated in SNRs. Calcula- 
tions performed within kinetic nonlinear theory demonstrate that the expected 
GCR spectrum, produced in SNRs, in satisfactory way fits the existing GCR 
data up to the energy 10^"^ eV. The first knee in the observed all particles 
GCR spectrum is attributed to the maximum energy of protons, produced 
in SNRs. The steepening of the all particle GCR spectrum above the knee 
energy 3 x 10^^ eV is a result of progressive depression of contribution of light 
CR nuclei with the energy increase. Such a scenario is confirmed by KAS- 
CADE experiment which shows r elatively sharp cutof f spectra of GCR species 
at energies emax ~ 3Z x 10^^ eV (lAntoni et all , boosi ). 



The nature of the whole GCR spectrum, which does not require any signifi- 
cant contribution of galactic CR sources other then SNRs, can be the following. 
GCR spectrum up to the energy e = 10^'' eV is dominated by the contribution 
of galactic sources, whereas at e > 10^^ eV GCRs are predominantly extra- 
galactic. Such a picture is consistent with the analys is, presented here, and 
with the analysis of ultra high energy CR properties (lAloisio et al.l . |2006| ) as 
well. Detail experimental study of GCR spectrum and composition could find 
whether indeed the transition from galactic to extragalactic component takes 
place in energy range 10^^ — 10^^ eV. 



Due to high interior magnetic field in young SNRs the 7r°-decay 7-rays gener- 
ated by the nuclear CR component as a rule dominates over 7-rays, generated 
by electron CR component, and calculated 7-ray fiux fits existing data. 



New generation of ground-based 7-ray detectors CANGAROO III, H.E.S.S., 
MAGIC, VERITAS will presumably measure significantly more accurate 7-ray 
spectra from already detected SNRs, hke Gas A, RXJ1713.7 and Tycho's SNR, 
as well as from a number of new sources. In this regard considerable interest 
as a potential 7-ray source s represent Tycho's SNR, SN 1987A and Kepler's 
SNR (IBerezhko et al.l . 120061 ) eventhough the expected 7-ray fiuxes from these 
objects only slightly above the sensitivity of modern TeV 7-ray instruments. 
Together with data, which soon will be obtained by space 7-ray observatory 
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GLAST, they will reproduce the full 7-ray spectra of some number of individ- 
ual SNRs, as well as the Galactic background spectrum. Although the existing 
data together with their theoretical interpretation presents a consistent pic- 
ture of GCR origin, these new measurements will allow for detailed study of 
SNRs and the CR production inside them. 
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Fig. 1. (a) Shock radius i?s) contact discontinuity radius i?c, shock speed Fs, and con- 
tact discontinuity speed 14, for Tycho's SNR as functions of time, inclu ding particle 



accel eration; (b) total shock (cr) and subshock (fis) compression ratios (jVolk et al 
20061 ^. The dotted vertical line marks the current epoch. The solid and dashed lines 
correspond to the internal magnetic field strength = 240 /iG and (i?d = 360 /iG), 
respectively. T he observed mean size and speed of the shock, as determined by radio 
measurements (jTan and Gulll . Il985l ) , are shown as well. 
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Fig. 2. The ratio Rc/Rs of the radii of the c ontact discontinu ity and the forward 
shock of Tycho's SNR as a function of time ( Volk et al. . 20061 ). Sohd and dashed 
hues correspond to the sa me two cases as in Fig.l. The experimental point is taken 
from lWarren et aD ( 2005 ). 
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Fig. 3. Spatially integrat ed spectral energy distribut ion of RX J1713. 7-394 6 



teerezhko and Volkl . l2006l^. The ATCA rad io data (cf. lAharonian et al.l. l2006al') 



ASC A X-ray data (cf. Aharonian et al. . 20051 ). and H.E.S.S. data ( Aharonian et al 



20071'! are showii . The EGRET upper limit for the RX J1713.7-3946 position 
( Aharonian et al. . 20051 ) is shown as well. The solid curve at energies above 10'' eV 
corresponds to 7r*^-decay 7-ray emission, whereas the dashed and dash-dotted curves 
indicate the inverse Compton (IC) and Nonthermal Bremsstrahlung (NB) emissions, 
respectively. 
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Fig. 4. Projected radial profile of the X-ray synchrotron em ission for the energy 
e^ = l keV from RX J1713.7-3946 (|Berezhko and Void , [iooi ^ . The solid line corre- 
sponds to the high-injection model; the dashed line again represents t he above pro- 



file, b ut smoothed to the resolution of the XMM-Newton data used in iHiraga et al 
the dash-dotted line corresponds to the test-particle limit. 




Fig. 5. The projected radia l profile of 7-ray emission for 7-ray energies = 1 TeV 
(|Berezhko and Volkl . boOfil ). The calculated radial profile is giyen by the solid line; 
the dashed line represents the calculated profile smoothed with a Gaussian point 
spread function of the width = 0.05°. For purposes of presentation both profiles 
are normalised to their peak values in this figure. 
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Fig. 6. (a) Shock radius i?s, shock speed Kj, gas density Ng and upstream magnetic 
field i?o at the current shock position; (b) total shock (a) and subshock (cts) com- 
pression ratios; (c) kinetic energy of ej ecta Ee,j and accelerated CR energ y content 
Ec as functions of time for SN 1987A (jBerezhko and Ksenofontovl. 120061) . Th e ob- 
served radius of the SN shock, as determined by radio (lManchesterelay,B and 
X-ray measurements (jPark et all 12004 ). are shown by circles and stars respectively. 
The scaling values are i?. = i?^ = 3.1 X 10^^ cm and Vi = 28000 km/s. The doited 
vertical line marks the current epoch. 
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Fig. 7. Sy nchrotron energy spectrum of SN 1 987A, calculated for the fiye eyolution- 
ary epochs iBerezhko and KsenofontoyI (|2006l ). The ATCA r adio dManchester et aU 



2002; Manchester and Gaensleil . I2OO1I ) and Chandra X-ray (jPark et alJ . I2OOJ) data 



for three epochs are shown as well. Higher measured fluxes correspond to later 
epoch. 
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Fig. 8. En ergy spectra for five group s of abundant GCR nu clei and all particle s 
spectrur n (iBerezhko and Volkl. [20071). The data of ATIC-2 ("Panov et all I2OO6I ). 
JACEE (jTakahashi et al.l . Il998l ) and KASCADE (jAntoni et al.. . .2005 ) experiments 
are shown as well. 
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